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A  THEORY  OF  THE  SHOCKWAVE 
PRODUCED  BY  A  POINT  EXPLOSION 


by 


Hans  G.  Snay 


ABSTRACT:  Approximate  integrals  of  the  spherical  blast 
equations  with  variable  lsentropic  exponents  are  derived. 
The  distributions  of  velocity,  density  and  pressure  within 
the  sphere  of  disturbance  are  expressed  in  polynomials. 
These  are  used  to  calculate  the  "reduced  energy"  vihich  is 
closely  related  to  the  relationship  between  shockwave  peak 
pressure  and  distance.  For  very  large  distances,  the  well 
known  asymptotic  behavior  of  weak  shockwaves  is  obtained. 
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This  report  gives  a  detailed  solution  of  the  problem  of  the 
shockwave  produced  by  a  point  explosion  in  water.  This  new 
approach  can  be  applied  to  other  media  and  is  generally  valid 
over  a  very  wide  range.  It  is  the  belief  of  the  author 
that  the  methods  described  represent  a  considerable  advance 
in  the  theoretical  treatment  of  blast  phenomena.  The  work 
was  done  under  Task  NO  701-267/76001/01040,  and  contributes  to 
the  solution  of  Key  Problem  No.  20,  Chapter  II,  "Key  Problems 
in  Explosives  Research  and  Development,  Part  II" ,  NAVORD 
Report  4299  (SECRET  Restricted  Data). 

The  author  wishes  to  thank  Mr.  J.  P.  Butler,  formerly  In  the 
Underwater  Explosions  Division,  for  his  help  in  the  preDaration 
of  this  report. 
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A  THEORY  OP  THE  SHOCKWAVE 
PRODUCED  BY  A  POINT  EXPLOSION 

CHAPTER  I  INTRODUCTION 

The  theoretical  treatment  of  shockwaves  caused  by  con¬ 
ventional  explosions  Is  complicated  by  the  presence  of  two 
media:  the  reaction  products  of  the  explosive  and  the  ambient 
medium  Into  which  the  shock  Is  transmitted.  Assumption  of  an 
Infinitely  small  "explosive  charge"  which  delivers  a  finite 
energy  removes  the  difficulty  of  the  two-medium  problem  and 
permits  study  of  the  shockwave  propagation  through  a  homogene¬ 
ous  medium.  Such  an  analysis  is  applicable  to  atomic  ex¬ 
plosions  or  to  electric  discharges  of  high  energy  In  small 
spark  gaps. 

Several  attempts  have  been  made  to  attack  this  problem 
[a]  -[e]**  the  best  known  is  Sir  Geoffrey  Taylor* s  now  classi¬ 
cal  treatise  "The  Formation  of  a  Blast  Wave  by  a  Very  Intense 
Explosion"  [a].  Recently  Lockwood  Taylor  [b]  has  shown  that 
the  solution  of  Taylor* s  differential  equations  can  be  given 
in  closed  form.  These  theories  apply  to  very  high  shock  pres¬ 
sures  only,  i.e.,  to  the  Initial  phase  of  the  blast. 

1 

The  present  paper  describes  an  approach  to  the  general 
problem  llch  also  holds  for  lower  pressures.  Particular  con¬ 
sideration  was  given  to  underwater  explosions.  Since  the 
thermodynamic  data  for  water  cannot  be  expressed  In  simple 
equations,  once  a  large  pressure  range  Is  considered,  the  cal¬ 
culations  were  set  up  for  the  use  of  tabulated  thermodynamic 
data.  The  methods  of  calculation  which  were  developed  can  be 
handled  with  ordinary  desk  computing  machines. 


#Such  letters  refer  to  the  list  of  references  at  the  end  of 
this  report. 
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CHAPTER  II  THE  TRANSFORMATION  OF  THE 
HYDRODYNAMIC  DIFFERENTIAL  EQUATIONS 

2.1  The  Hydrodynamic  Equations.  Even  for  strong  energy  dis¬ 
charges,  radiation  effects  are  unimportant  under  water,  because 
the  mean  free  paths  of  photons  and  electrons  are  small  In  media 
of  high  density.  Also,  it  is  well  known  that  viscosity  and 
heat  conduction  do  not  affect  the  fluid  motion  we  are  consider¬ 
ing  here  except  at  shock  fronts.  There,  these  effects  are  im¬ 
plicitly  accounted  for  by  the  Rankine-Hugoniot  shock  relations. 
Thus,  for  regimes  behind  shock  fronts  the  energy  is  transmitted 
by  pressure  forces  only  and  each  particle  has  constant  entropy 
as  long  as  it  is  not  traversed  by  a  further  shock.  This  situ¬ 
ation  is  described  by  the  well-known  hydrodynamic  equations. 


For  the  case  of 

spherical 

symmetry. 

these  are 

• 

• 

(2.1) 

u 

+  p  fir 

-  °  / 

(2.2) 

* 

«  pr  + 

J  y- 

O  . 

/» 

(2.3) 

S*  * 

O . 

Introducing  the 

isentropic  exponent 

(2.4) 

/  - 

/fAAJ  = 

( ?  -£jr/s 

j 

(2.2)  and  (2.3) 

can  be  combined  as  follows : 

(2.2a)  ~  prj  -  /  (fi±  +  +  ty) . 

The  symbols  are  defined  in  the  table  at  the  end  of  this 
paper.  It  should  be  noted  that p  designates  the  excess  pressure 
above  the  static  pressure  ^  .  This  convention  affects^  which 
is  here  defined  somewhat  differently  from  the  common  isentropic 
exponent,  ^  .  The  latter  is  related  to  £  as  follows: 
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(2.5) 


At  high  pressures,  there  is  no  difference  between  these  two 
magnitudes  but  for  p+o  ,  »,  whereas  remains  finite. 

2.2  The  Ranklne-Hugonlot  Conditions.  For  a  shock  advancing 
into  an  undisturbed  medium,  the  following  conditions  must  be 
observed : 


(2.6) 

(2.7) 


(2.8) 


-  fit  1 

u  = 


M, 

-  7%  ' 


E,  -E„  *  A  t  fit  -j. 


(The  subscript  1  designates  the  state  just  behind  the  shock 
front;  zero  refers  to  the  undisturbed  state). 

2.3  Thermodynamic  Data.  Equation  (2.8)  must  be  evaluated  using 
the  thermodynamic  data  of  the  fluid  considered.  The  evaluation 
is  simplified  by  the  introduction  of  the  Mreduced  internal 
energy” 


(2.9) 


J  -  3J 


EJs  p 

fi  P 


which  will  also  be  used  in  other  connections  in  this  paper. 

J  is  a  dimensionless  magnitude  related  to  the  heat  capacity. 
(For  an  ideal  gas  at  high  temperature  and  pressure, J  •  Cr/R 
where  cv  is  the  heat  capacity  at  constant  volume  and  R  is 
the  gas  constant).  The  relation  between  J  and  p  is: 


(2.io)  ^  -  -j-  ■ 
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This  relation  holds  generally  for  any  type  of  fluid.  For 
constant  heat  capacity,  the  last  term  vanishes  and  one  obtains 
a  relation  which  is  widely  used  in  the  literature. 

Once  Jffyp)  is  known,  one  can  readily  determine  the  shock 
pressure  as  a  function  of  the  density  p,  and  the  magni¬ 
tude 

<2-4a) 

This  is  analogous  to  the  isentropic  exponent  (2.4)  except  that 
the  differential  quotient  is  taken  along  the  Rankine-Hugoniot 
adiabatic  instead  of  the  isentropic.  According  to  the  Rankine- 
Hugoniot  conditions,  ' /  fyn  vanishes  for  infinitely  high 

pressures  whereas  remains  finite. 

2.4  Change  from  Free  Boundary  Conditions  to  Fixed  Boundary 

Conditions 

Equations  (2.1)  to  (2.3)  constitute  a  free  boundary  problem 
since  the  boundary  conditions  must  be  fulfilled  along  an  un¬ 
known  line,  -  namely,  the  radius-time  curve  of  the  sphere 
bounded  by  the  shock  front.  The  hydrodynamic  equations  will 
now  be  transformed  in  such  a  way  that  three  partial  differential 
equations  are  obtained  whose  boundary  conditions  must  be  ful¬ 
filled  at  a  fixed  and  known  place.  This  is  possible  by  adding 
one  more  equation,  namely  a  relationship  for  the  position  of 
the  shock  front. 

The  transformation  is  readily  performed  by  the  introduction 


of  the 

following 

reduced  variables 

• 

• 

(2.11) 

Reduced 

velocity: 

(2.12) 

Reduced 

density: 

x  - 

f, 
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(2.13) 

Reduced  pressure:  ^  m 

(2.14) 

Reduced  distance:  J  * 

K 

All  magnitudes  with  the  subscript  1  refer  to  the  shock  front 
and  are  functions  of  time  only;  in  particular,  we  note  that 

<dLr,/dt  -  U  . 

With  the  use  of  these  variables  and  the  Ranklne-Hugonlot 
relations  (2.6)  and  (2.7)#  the  hydrodynamic  equations  take  the 
form: 

(2.15) 

r*  *  fa 

+  (#* 

- 

(2.16) 

(2.17) 

y  Y  ft  t 

A- A 

Y  “  f  /('fat  ¥) 

+  m)*  , 

+  zftt  ( to>  V* 

j  pfd* 

m 

We  Introduce  the  abbreviations: 

(2.18)  H,  -  £  dfe*L> 

• 

(2.19) 

JjL  /4C  / 

^  ~~s* 

-  fc*  tfaj' 

(2.20) 

G,  =  JL  p, 

**• 

-  H>(  fan  • 

Since  fa*  is  known  as  a  function  of  />,  ,  only  H,  remains 

to  be  determined.  (This  will  be  done  in  the  next  chapter). 
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Using  &, ,  we  can  eliminate  the  time  £  from  the  partial 
differential  equations  and  obtain  equations  with  and  fi,  as 
Independent  variables ,  where  p,  takes  the  place  of  £  ; 


(2.21) 


(2.22) 


(2.23) 


^  - 


%  = 


JL  L>  * 

fa  +  -j-  6i  Q_  +  P'  b  ^ 

JL  MP*  +  if)  t—th—L  ^adJjjdp. . 

I  -JL  -  JL 

A-fi.  f 


JL 

j 


_  M- 

T 


-  SL 

3 


2.5  Boundary  Conditions.  The  transformation  of  the  partial 
differential  equations  has  simplified  the  boundary  conditions. 
These  are: 

<f  -  /  ) 

(2.24)  X  *  I  i  /«•  I* 

f  ~  I  ) 

At  the  center  of  the  sphere,  J  -  0,  the  velocity  must  vanish. 
Hence , 


(2.25)  f  -  0  for  f  m  0. 

Using  (2.1),  it  can  be  deduced  that 

(2.26)  ^  -  °  for  ^  -  0^ 
if  X  a°d  ^ 

The  partial  differential  equations  (2.21)  through  (2.23) 
have  some  unique  properties  when  combined  with  the  boundary 


It  are  finite  or  zero  at 

J  3 


0. 
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conditions  (2.24).  The  derivatives  with  respect  to  p,  vanish 
at  »  1.  This  means  that  the  derivatives  of  f  >  X  an<*  ¥* 
respect  to  £  at  £  -  1  are  functions  of  fl,  and  Ht  only, 
provided  the  thermodynamic  properties  of  the  medium,  expressed 
by  /  and  fTtH  ,  are  given.  This  also  holds  for  the  higher 
derivatives.  The  consequences  of  this  behavior  will  be  dis¬ 
cussed  in  the  next  chapter  after  ths  relation  fcr  H,  has  been 
found  • 


CHAPTER  III  GROSS  PROPERTIES  OP 
THE  SPHERE  OP  DISTURBANCE 

The  shockwave  which  is  emitted  from  the  point  explosion  is 
spherically  symmetric.  We  call  the  region  which  is  enclosed  by 
the  shock  front  the  "sphere  of  -dis turbance" . 

3.1  Average  Density  Relation.  If  we  neglect  the  presence  of 

the  device  which  discharges  the  energy,  the  medium  inside  the 

sphere  of  disturbance  is  the  same  as  that  outside  the  sphere. 

Since  the  ambient  medium  remains  undisturbed  until  the  shock 

front  arrives,  the  average  density  within  the  sphere  must  be 

that  of  the  undisturbed  medium,  or 

/ 

(3.D  -  {*> 

o 

This  relation  is  consistent  with  the  hydrodynamic  equation  for 
the  conservation  of  mass.  It  makes  a  simple  statement  about 
the  behavior  of  X  •  Similar  relations  for  (p  and  ^  do  not 
exls  t . 

3.2  The  Energy  Integral.  In  order  to  derive  an  expression  for 
Hi  ,  we  consider  the  total  energy  of  the  sphere  bounded  by  the 

shock  front.  As  this  sphere  expands,  its  energy  in  excess  of 
that  of  the  undisturbed  medium  must  remain  constant  and  equal 
to  the  total  energy  Q  of  the  initial  discharge. 
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The  total  energy  Increment  of  the  sphere  la 


r, 


K 


/t  2 

(J&i  +  p£)rldr  -  in-  //4  £.  r  dr, 


where  the  first  term  Is  the  sum  of  the  kinetic  and  Internal 
energies  inside  the  sphere  and  the  last  term  is  the  total 
energy  of  the  undisturbed  fluid  contained  in  a  sphere  of  equal 
size. 

Introducing  the  dimensionless  variables  p  ,  X  »  P  and  the 
magnitude  J  ,  (2.9) ,  we  obtain 


(3.3) 


i 

Uj<xr‘ *  rjfif 

O 


The  last  term  vanishes  due  to  (3.1) .  The  first  Integral  is 
abbreviated,  using  (2.6),  by 

/ 

(3.4)  t,  -  Jjfxr1  - 


This  magnitude  represents  the  reduced  energy  of  the  sphere  of 
disturbance.  The  total  energy  is 


(3.5) 


Equation  (3.3)  provides  a  convenient  relationship  between  the 
shock  pressure  and  the  distance  A  for  any  given  energy  yield 
Q  .  Hence,  the  primary  task  is  to  find  as  a  function  of 

p>  • 
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Setting  a?  ■  »'  “  obtain  after  rearrangement 


H, 


(3.6) 


'  +  fa  l#*  + 

A 


-3 


<4  -/V 


'  *  +  ftl 


A 


Z/T3£a> 

the  aesl] 


where 

(3.6)  together  with  (3.4)  gives  the  desired  relationship  be¬ 
tween  /V,  and  p4  .  The  solution  of  our  problem  must  satisfy 
the  three  partial  differential  equations  (2.21)  through  (2.23) 
as  well  as  the  integro-differential  equations  (3.6)  and  (3.4). 

It  is  now  possible  to  see  the  significance  of  the  trans¬ 
formation  (2.15)  to  (2.17).  We  assume  that  p  and  are 
known.  As  shown  in  the  previous  chapter,  the  first  and 
higher  derivatives  of  ff,  y  and  ^  with  respect  to  J  ,  at 
-  1,  functions  of^  and  J$  only.  If  if9  J(  and  ^  were 
analytic  functions,  they  could  be  expanded  around  f  *  1,  in 
series.  Introducing  these  into  the  Integral  for  £  ,  (3.4), 
and  using  (3.6),  we  obtain  an  ordinary,  first  order  differ¬ 
ential  equation  for  .  Thus,  our  transformation  and  the 
expansion  around  £  -  1  reduces  our  problem  to  one  ordinary 
differential  equation.  Once  {//>,)»  or »  is  known,  the 
problem  is  solved. 

In  this  paper  we  will  not  discuss  whether  or  not  this 
type  of  solution  is  possible.  For  practical  calculations, 
expansions  of  this  kind  are  usually  very  inconvenient  if  a 
large  number  of  terms  is  involved.  However,  there  remains 
the  highly  significant  point  that  the  behavior  of  our  three 
reduced  variables  near  J  «  1  gives  the  greatest  contribution 
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to  if,  ,  (due  to  the  factor  J2  which  occurs  under  the  integral) 
and  that  this  behavior  can  be  calculated  without  solving  the 
partial  differential  equations. 

Closely  related  to  //,  la  a  magnitude  of  great  practical 
Interest  In  any  shockwave  problem,  namely,  the  inclination  of 
the  pressure-distance  curve  In  a  logarithmic  A  *  ^  plot: 


(3.7)  -  JL  U  -  -  -r— — - 

This  shows  that  for  ■  0,  the  peak  pressure  is  Inversely 

proportional  to  the  third  power  of  the  distance.  This  holds 
for  very  high  pressures.  For  large  distances  £  and  low 
pressures,  the  wave  assumes  a  nearly  acoustic  behavior  and 
approaches  unity,  as  will  be  shown  In  Chapter  VII. 

For  a  finite  fa*  and  negative  ft  ,  the  pressure  may  decrease 
even  more  rapidly  than  with  the  cube  of  the  distance. 

3*3  The  Dissipated  Energy.  At  the  shock  front  mechanical 
energy  is  constantly  dissipated  and  converted  into  irreversi¬ 
ble  energy.  According  to  the  second  law  of  thermodynamics, 
the  Increase  of  dissipated  energy  Is 


(3.8) 


4nrfa  f  r*(S)d.S 

Sfo) 

4  j rr,  lp0  h 


T(S)  is  the  lowest  permanent  temperature  which  can  occur  In 
the  system.  For  liquids,  this  temperature  corresponds  closely 
to  that  which  the  fluid  attains  after  passage  of  the  shock- 
wave  and  return  to  zero  excess  pressure. 

The  Integral  In  (3*8)  Is  commonly  called  the  dissipated 
enthalpy  Increment  h  •  If  the  lsentropic  p-p  relationship  Is 
known,  h  can  be  conveniently  determined  from  the  difference 
of  the  enthalpy  Increments  at  the  shock  front,  ,  and  that 
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of  an  lsentropic  expansion  to  zero  excess  pressure: 

A 


(3.9) 


A  - 


const ) 


AH  -  (~fr.  +  -fsr,  )  ■ 

This  relation  shows  that  h  Is  a  shock  parameter  like  U  or 
M, ,  and,  therefore,  la  a  function  of  fi,  only.  Behind  shock 
fronts,  h  Is  constant  along  any  lsentropic  line. 

The  amount  of  dissipated  energy  which  is  within  the  sphere 
of  disturbance  can  be  found  by  integrating  (3*8)  along  the 
shock  front.  Introducing  defined  by 


(3.10)  Q-D'm  -  ^  r‘3 & 

the  following  expression  is  found: 


a  9 


(3.11) 


soL£*** 


=-j- 


)  - 


Zp±  hfp,) 

7z*  i  P>  (P>  -?•) 


J+  folfa.  ~ 


This  and  (3-7)  establish  a  relationship  between  £  and  2d*sj  • 
Either  one  of  these  equations  can  be  used  to  determine  /J,  . 

It  is  convenient  to  use  ft  for  high  values  of  p,  ,  whereas, 
(3.11)  will  be  used  for  the  case  where  pj  approaches  ,  i.e. 
for  very  low  shock  wave  pressures. 


CHAPTER  IV  THREE  IMPLICIT  INTEGRALS 


In  this  chapter  we  will  derive  implicit  integrals  of  the 
three  partial  differential  equations  (2.21)  through  (2.23). 
These  integrals  contain  derivatives  with  respect  to  p,  and 
are  therefore  not  complete  solutions  of  our  problem.  But  they 
give  the  initial  distribution  of  our  solution  and  provide  an 
understanding  of  the  behavior  of  the  general  solution. 
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Further,  they  can  be  used  for  iterations  to  improve  approximate 
solutions. 

4.1  Integral  Obtained  From  Energy  Considerations.  Consider  a 
sphere  lying  concentrically  within  the  sphere  of  disturbance. 

The  inner  sphere  does  the 

t 

work  ■  4/r  J  ^ 

t, 

on  the  surrounding  medium  while  it  expands.  It  contains  the 

r' 

energy-  [(•%?+ 

/ 

In  these  relations,  a7  denotes  the  radius  of  the  inner  sphere. 

According  to  the  assumption  made  in  paragraph  2.1,  the  change 

of  energy  of  the  sphere  is  equal  to  the  rate  of  mechanical 

work  done  on  the  surrounding  medium: 

r‘ 

(4.1a)  -  f/j*1  +  EpJrUr  -  + 

o 

or,  since  sdLr'J/dbt  » 

r  * 

(4.1b)  JipjptoE  +  p(£-fe)Jr*J  dr  +  r  + p(£-Fm)  *  -  0 

0 

In  (4,lb)  E  has  been  replaced  by  E-EQ  since,  by  virtue  of 

(2.2) 


The  partial  derivative  in  the  integral  of  (4.1b)  refers  to 
constant  r  .  If  we  change  to  the  variable  we  must  note  that 


~  n't 

@f)*  ■ 
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Introduction  of  the  reduced  variables  yields 

S' 


(*.2) 


Introducing  further  J/ 

(4.3)  ?  -  -  3  for*  + 

o 

and  making  use  of  (2.19) >  (2.20)  and  (3.6),  we  obtain: 


=  0 


firat  term  of  (4.2)  -  ^  £  ?£  +  £■  ^ 


(4.4) 


■  -  J£1 


■  f  M**  - 

Rearrangement  of  (4.2)  and  (4.4)  glvea  finally 


(4.5)  » 

X<fl 


A-/>. 


r/y 


W  [('*  t  '  fgf-zk]- 


In  this  equation  the  prime  Is  omitted  on since  It  now  repre¬ 
sents  the  Independent  variable, 

4.2  Integral  for  the  Reduced  Velocity.  By  combination  of  (2.21) 
and  (2,23)  we  obtain 


•  l.  **a#*jSFj^ft»***‘«*] 


(4.6) 
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Prom  (4.5), 


(4.7) 


JL  gg 

^  7&.  ~  ** 


»/? _ 


Jf  -9  ' 

where  we  have  used  the  following  abbreviations: 


(4.8) 


$  - 

?  -  ta*>) 

o  -  ja  j  -  r 

r  '  JL  Jk  '  fr/7 


Furthermore,  we  abbreviate 


3^ 


(4.9) 


H  - 


+  6.  72=?, 
M  *  6- 

&  +  6-  teS, 


ro 


(4.6 )  then  takes  the  form: 


(4.10) 


=  (3j+})pl  +(H  +  ~  dis¬ 


integration  gives: 


(4.11) 


,«  f  (±ll!l  tJtlL  4± 

"“J  '  J  ('V+?7^*  +  r«  «•**-  7 


co»s+. 


If  the  coefficients  of  ^  in  (4.11)  are  constants,  the  integral 
can  be  expressed  in  closed  form: 

14 
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if  .  &!•*  +  j  **  J^(  Ml  zJ± ) 

xrmrj (  /- 1 1 


(4.12) 

i  KoSi  4*  KSl  t  Kx  £ 

ii  Cit-S. ) 

where 

(4.13) 

K,  -  (ft  +  D/Pw) 
Kt  =  l(4r*p) 

and  where  S, ,  &x  are  the  roots  of  the  quadratic  equation  in  the 
denominator  of  the  integrand. 

4.3  An  Integral  Involving  the  Reduced  Density  and  Pressure, 
Prom  the  differential  equations  (2.22)  and  (2.23)  one  obtains, 
using  the  abbreviation  (4.9)f 


(4.14) 


-  -  a,  %  , 


Integration  ylelda 


+  - All 


.  A  -?f  } 

*  t 

It  is  interesting  to  note  that  both  integrals  remain  finite  for 
all  values  of^  from  zero  to  unity. 

For  '/fan  m  o  and  constant^  the  expressions  (4.5),  (4.12) 
and  (4.15)  correspond  to  the  Lockwood  Taylor  solutions  of  the 
strong  point  blast  wave  [b]. 


15 
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CHAPTER  V 
THE  TAYLOR  CASE 

5.1  The  Initial  Distribution  of  the  General  Solution.  0.  I. 
Taylor  has  described  a  special  case  which  applies  to  the  early 
phase  where  the  shock  pressure  is  very  high  [a].  Corresponding 
solutions  of  the  hydrodynamic  equations  are  called  nprogres8lve 
waves'1  by  Courant  and  Friedrichs [f] .  The  essential  feature  of 
this  solution  is  that  the  velocity ,  pressure  and  density  distri¬ 
butions  inside  the  sphere  do  not  Change  shape  as  the  shockwave 
advances;  only  the  scale  changes.  In  our  coordinates,  this 
means  that  the  p-jf,  Jf-jf  and  y'-J -curves  remain  the  same;  they  do 
not  change  with  time  nor  with  ft  which  replaces  time  in  our 
analysis . 

A  glance  at  equations  (2.21)  to  (2.23)  shows  immediately 
that  this  can  be  true  only  if  ft  is  constant.  Constancy  of  />, 
is  also  required  in  Taylor's  treatment  and  is  Justified  for  high 
pressures  where,  according  to  the  Ranklne-Hugonlot  conditions, 

becomes  independent  of  ftt  i.e.,  '/ft h  *  0.  This  reduces  the 
partial  differential  equations  to  ordinary  differential  e- 
quations . 

Since  '/ft*  -  0  for  infinitely  high  pressure,  the  Taylor 
solution  provides  the  initial  distribution  from  which  the 
general  solution  of  the  partial  differential  equations  develops. 
Thus,  we  assume  that  we  have  Taylor  distributions  from  the  very 
beginning  and  that  these  remain  stationary  as  long  as  does 
not  change. 

(Actual  explosions  do  not  have  such  distributions  from  the 
very  beginning.  However,  it  is  safe  to  assume  that  the  distri¬ 
butions  ef  pressure,  velocity  and  density  will  quickly  converge 
to  the  Taylor  distributions  or  the  subsequent  distributions 
calculated  on  the  basis  of  an  initial  Taylor  solution.) 


16 
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5.2  The  Taylor  Solution,  Por  a  true  point  explosion  both  the 
pressure  and  the  temperature  are  infinitely  high  initially, 
whereas  the  density  is  finite.  Under  such  circumstances  the 
medium  will  be  completely  ionized,  l.e.,  all  the  atoms  of  the 
medium  will  be  completely  stripped  of  their  orbital  electrons. 
Even  at  relatively  high  densities  such  a  medium  behaves  like  an 
ideal  monatomic  gas.  Effects  not  present  in  the  ideal  gas  theo¬ 
ry  are  those  associated  with  the  electrostatic  forces  between 
atoms  and  the  radiation  pressure.  We  neglect  both  of  these  in 
this  analysis.  Then  <f  is  constant  throughout  the  sphere  and 
equal  to  5/3 *  the  value  for  an  ideal  monatomic  gas.  Further¬ 
more,  we  have  the  following  well  known  relations  which  hold  for 
ideal  gases  with  constant  heat  capacities: 


(5.1) 


/.  ■/  J 


With  f  ■  5/3#  we  have  p/  -  ^po  .  This  is  the  highest  value  of 
P,  which  can  occur  in  our  problem.  Por  the  Taylor  case,  we 


have 


(5.2) 


JL 


0 


The  Integral  (4.5)  becomes: 
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JL 

(5.3) 


A  zJS 
r*  -  fa. 


The  equation  for  the  reduced  velocity  ia  also  simplified  since 
the  coefficients  of^  in  (4.11)  are  constants.  With 

(5.*)  ?  -  -  -£#> 

and  j  *  p,Jpa  , 

we  obtain  from  (4.13) • 


(5.5) 


Ko 

K, 


_ £l- 

J-JjL.  Jr  _BL 

+  t  P.~f» 


K 


The  roots  of  the  quadratic  are 

s,  . 

'  £a£js, 


(5.6) 


fi 


If  we  examine  (2.23)  at 


(5.7) 


0,  we  find,  since  vanishes. 


We  will  frequently  deal  with  this  magnitude  and  use  the 
abbreviation 


(5.8) 

The  equation  for  the  reduced  velocity  then  becomes: 
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m 


(5.9) 


^  _  j, 

j  - c*  '  -  </v 
-  &  +*)-' 


i£t 


-L  A^Je. 

S 


3/t 

9  =  -  ffi,  jfefe,  s' P' 

df,  +  Pt> 


Finally,  we  obtain,  from  (4.15),  t.he  simple  expression 


(5.10) 


Combination  of  (5.3)  and  (5.10)  permits  the  computation  of 
J'  and  J(  ,  once^/jf'  has  been  found  from  (5.9).  All  three 
functions  are  equal  to  unity  forj?  -1.  p  and^  vanish  at 
f  "  0,  whereas^  remains  finite.  It  is  of  interest  to  study 


the  behavior  of  these  three  functions  at  f  »  0. 

Prom  (5.9),  it 

can  be  deduced  that 

(5.11) 

jf-  o 

~^*)]  ~  .?  i 

where 

z 

(5.12) 

v  [  cy  •  i:  -]  'c- 1 

Hence,  we 
of  f  is 

see  that. 

for  0  in  the  Taylor  case,  the  behavior 

(5.13) 

9  - 

c°f  r  c'$  * . 

or,  for/ 

-  5/3, 

(5.13a) 

9  = 

„  7r 

#<Pf  <?<?<??  £  +  . 
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Combining  (5.3)  and  (5.10),  we  obtain  for  the  reduced  preasure 

V  » 


If 

(5.14) 


This  shows  that  %  is  finite. 
Putting  (5.13)  inco  (2.23),  we 

(5.15)  y  =  Ao  + 


(%  =0.3060  for  r  -  5/3). 


obtain  for  small  values  of  £ 


where 

(5.16) 


Hence,  for  ■  5/3,  we  have  near  ^  *  0 


s:r 


03060  +  + 


(5.17)  f  ~ 

For  the  reduced  density,  we  obtain  directly  from  (4.15): 


3 


(5.18)  =  s.i 


with 

„  4  ^  AA- 

-  i3'  "A  <T0 " 

-X 

=  •  2  A 

/*•  J 

1 

(5.19) 

^  <?.  J* 

A- 

0  • 

5.3  The  Primakov  Case.  For  lf  the  Taylor  solution  takes 
on  an  Interesting  form,  first  found  by  Primakov  [g].  In  this 
case,  p,  »  ^  P*  and  C0  ■  1.  C,  ,  >4,,  A,  and  3* vanish .  The 
solution  of  (2.21)  to  (2.23)  is 


(5.20) 


Jr  *  /  A-  /'*  /  • 
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Of  course,  this  solution  fulfills  the  average  density  con¬ 
dition  (3.1) •  The  Primakov  case  Is  often  mentioned  in  con¬ 
nection  with  point  explosions  under  water,  because  an 
lsentroplo  pressure-density  relationship  using  the  exponent  7 
gives  a  good  representation  of  the  behavior  of  water.  This, 
however,  holds  for  moderately  high  pressures  only.  This  Is 
evident  If  we  consider  the  highest  density,  which  Is  1.333 /?© 
In  the  Primakov  case  but  4^#for  a  completely  dissociated  and 
Ionized  medium. 


CHAP1ER  VI  METHOD  OP 
ATTACK  FOR  THE  GENERAL  CASE 

6.1  The  Determination  of  In  order  to  perform  calculations 
using  the  transformed  differential  equations  (2.21)  to  (2.23), 
it  is  necessary  that  H,  be  known  as  a  function  of/^  . 

The  proposed  method  assumes,  as  an  Initial  step,  several 
arbitrary  values  for  ,  or  /3  ,  such  as  1,  0,  -1,  -2. 

Then  ,  L ,  and  can  be  computed  as  functions  of  p,  ,  since 

p, )  is  known.  Next,  ,J[  andp  must  be  determined  from 

the  partial  differential  equations  by  one  of  the  methods 
described  below.  Subsequently^,  Is  calculated  using  (3.4). 
When^  Is  plotted  versus^  —  P°)J P*  In  logarithmic  scale,  thb 
Inclination  of  the  derived  *},  curve  must  coincide  with  the 
value  of ^  which  was  used  to  calculate  this  particular  curve. 
Such  a  curve  is  readily  drawn  In  the  same  way  as  Is  the 
graphical  solution  of  differential  equations  by  the  method  of 
Isoclines.  Figure  1  Illustrates  the  procedure.  Usually,  the 
first  step  gives  satisfactory  accuracy.  The  procedure  can  be 
repeated  assuming  more  closely  spaced  values  ot  /3  ,  until  the 
desired  accuracy  Is  obtained. 

This  method  Is  practical  only  If  ^  can  be  determined  from 
computations  which  are  not  too  cumbersome.  This  Is  a  formi¬ 
dable  task,  if  exact  solutions  are  desired,  especially  for 
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♦ 


variable  f  and  7  .  However,  it  la  possible  to  devise  approxi¬ 
mations  which  simplify  the  computations  considerably  without 
too  great  loss  of  accuracy. 

6.2  Differential  Quotients  tor}  ■  1.  The  reduced  functions 
p  ,  X  and  p  have  the  value  one  at  jr  *  1.  Henoe,  the  deri¬ 
vations  with  respect  to  pt  vanish  In  (2.21)  to  (2.23),  and  we 
obtain  the  following  simple  relations: 


(6.1) 


(6.2) 

(6.3) 


The  second  derivatives  are: 


W,  v~  L,  * 

hzh  ( ?/+  Z  +  &,)  > 
[*(*>/+*)  *  ".] 


(6.4) 


(6.5) 


(6.6) 
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These  and  the  higher  order  derivatives  are  completely  determined 
as  functions  of  ft,  ,  once  H,  or  jS  and  the  thermodynamic  data  of 
the  medium  are  given.  The  case/V/ty^-^must  be  excluded. 

This  situation  is  approached  torft-+ft0  and  is  discussed  later 
in  Chapter  VII. 

6.3  The  Behavior  Kear  the  Center  of  the  Disturbance.  The 
Taylor  solution  shows  that  the  density  is  zero  at  the  center 
of  the  explosion  but  the  pressure  is  finite.  If  we  consider 
(2.9a)  for  the  center  (*c  ■  0),  we  have 


(6.7) 


For  a  constant  >  ■  x 


(6.8) 


P 

r 


this  becomes  the  familiar  lsentroplc 

_  <K 

*  c  P 

•  t’mzJ 

^  c  fi  r. 


The  lsentroplc  constants,  C  and  C.  \  are  infinitely  large  at  the 
center.  This  is  a  consequence  of  the  assumption  of  a  point 
explosion,  where  at  the  beginning  the  pressure  and  temperature 
are  infinitely  high;  it  holds  not  only  for  the  Taylor  solution 
but  also  for  the  subsequent  stages.  Therefore,  the  density 
remains  zero  and  the  temperature  infinitely  high  at  the  center 
as  long  as  the  basic  assumptions  of  our  calculations,  in  par- 
tlcular  the  absence  of  radiative  energy  transport,  are  valid. 
These  are  clearly  not  fulfilled  at  the  center  since  high 
temperature  and  low  density  necessarily  swan  strong  radiation. 
For  underwater  explosions,  this  inconsistency  with  the  physical 
picture  is  restricted  to  the  almost  empty  spaces  adjoining  the 
center.  At  some  distance  from  the  center,  radiation  effects 
will  become  small  because  of  the  Increased  density.  Actually, 
the  hot  center  will  cool  down  very  rapidly  due  to  radiation, 
equalizing  the  temperature  in  the  core.  This  temperature  may 
be  too  low  for  appreciable  radiation  to  the  surrounding  medium, 
but  still  high  enough  to  cause  dissociation.  This  may  Indicate 
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that  our  basic  equations  are  acceptable  approximations  and  that 
j9  ■  5/3  •  This  is  the  same  value  as  for  the  initial  distri¬ 
bution;  thus,  near  the  center^  is  constant  with  pt  or  with 
time. 

Prom  (4.15)  it  can  be  conclude'*  that 


(6.9) 


with 

(6.10) 


Inserting  (6.9)  into  the  partial  differential  equations  (2.21) 
to  (2.23)  one  finds  that  the  same  expansions  as  in  the  Taylor 
case,  (5*13)  and  (5*15),  satisfy  these  equations  near  $  »  0, 
namely 


(6.11) 

(6.12) 

(6.13) 


+  c'$  *•••• 


x 

f 


tr>  +  i. 

4o  +  + 


+  + 

-*•1  6.. 


This  expansion  breaks  down  if  Ac  becomes  zero  as  illustrated 
by  the  Primakov  case.  Therefore,  it  is  of  Interest  to  examine 
the  conditions  for  which  A,  may  vanish.  We  set 


(6.14) 


-Avyv  lL 

3 r*o  * 


Then,  according  to  (4.15) 


(6.15)  $+0 

Putting  this  Into  (2.21), 


* 
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shows  that  for/#  <  2,  A0  «  o  can  only  occur  if  we  permit  a 
singularity  in  ^at  the  center.  This  may  occur  when  standing 
waves  within  the  sphere  are  reflected  at  the  center  but  can  be 
excluded  for  any  other  condition.  Thus  the  expansions  (6.11) 
to  (6.13)  hold  for/;  <2. 

The  coefficients  A  ,  8  ,  C  in  (6.11)  to  (6.14)  are  functions 
of  fi,  and  are  determined  by  ordinary  first  order  differential 
equations.  The  first  four  of  these  are: 


(6.17) 

3/0C0  +  H,  +6' 

Am  /■> 

(6.18) 

-C.J 

(6.19) 

*•(&.  -  c‘)  -  *<■ 

-  6,  s.  it*  -  o 

(6.20) 

_ *  2  A  y 

- c-  *  s'  - 

o 


These  are  four  equations  in  five  unknowns ,  namely A*9At,  , 

C9  ,  Ct  .  Adding  one  more  differential  equation  always  brings 
in  one  more  unknown.  Hence,  this  expansion  around  J  ■  0  remains 
undetermined,  as  it  was  in  the  Taylor  case.  An  additional  ex¬ 
pression  which  relates  the  expansion  around  f  -  0  with  the 
complete  solutions,  is  required. 
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Several  such  expressions  can  be  set  up.  The  most  con¬ 
venient  one,  equation  (6.25) »  will  be  discussed  in  the  next 
paragraph.  Another  interesting  relationship  can  be  found  by 
combination  of  (6.17)  and  (6.19)  and  by  integration  from  the 
initial  density  to p(  : 

A  ,  */* 

(6.2i)  -  A -A  t.P  fyv  ) 

A.  (f>u)  P>  -A  %  ffr)  (  a  Z'ffi')/ 

This  Integral  requires  the  knowledge  of  %  and,  therefore ,  is 
less  convenient  to  use  than  equation  (6.25)* 

6.4  Polynomials  for  the  Reduced  Velocity,  Density  and 
Pressure.  Since  the  behavior  of  the  solutions  near  f  ■  1  and 
-  0  is  known,  one  may  construct  approximate  distributions 
f or  f ,  X  and^  by  merging  these  solutions  in  the  intermediate 
range.  One  of  the  simplest  ways  to  do  this  is  to  use  poly¬ 
nomials  in  J  ,  with  exponents  in  terms  of  at?  t  which  satisfy 
the  partial  differential  equation  near  »  0.  The  coefficients 
are  determined  in  such  a  way  that  the  polynomials  satisfy  the 
boundary  conditions  and  have,  up  to  a  certain  order,  the  correct 
behavior  near  f  «  1  and  f  •  0. 

If  only  first  order  terms  are  used,  we  have  four- term 
polynomials  for  j(  and  if  and  a  three-term  polynomial  for  : 


(6.22) 

r  - 

i  --  O  ° 

(6.23) 

x  = 

/:  e 

(6.24) 

f  - 

3  2  ) 

z.  i 

JS 
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-4  .  ' 


There  are  eleven  coefficients  to  be  determined  for  which  we 
have  the  following  conditions : 

(<?,  - 1 

(a)  Three  boundary  conditions  at  j-  •  1  -  X,  «  1 

U,-i 


(b)  Three  differential  quotients  at  j?  »  1 


i!1  (6.1) 
X,'  (6.2) 

V  (6.3) 


(c)  The  average  density  condition  (3.1) 

(d)  Three  ordinary  differential  equations  involving  4*  ,  Ba  , 

,  A  :  (6.17),  (6.19),  (6.20) 

(e)  One  expression  which  relates  the  expansion  around  £  *  0 
to  the  entire  solution,  preferably  a  condition  for  A*  . 

The  last  condition  (e)  can  be  found  from  (2.21).  Solving 
for  ^  yields 

(6.25)  P»  -  ,°>fSX  f  ~  X  ~ 

The  Integration  from  £  *  0  to  J  *  1  can  be  readily  performed 
with  the  use  of  the  polynomials  (6.22)  to  (6.24). 

With 


(6.26) 


we  obtain 


Sj  =  1  J± 

rJ  Uo  j  *  ‘ 


(6.27) 


%  “  A,  -  /-  C.(j3,  -  2- fix  *fi,)  -  Vj 

+  fa  { /K*  (A  -  A  j  -  j 

4  ^  fXCnCipi  +  (7, SC,  i-  /SCeCj  )pj 

+  2l.SC,Cifa  +  "Mr/  *  4?^ -J2& 
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(6.27)  is  an  ordinary  differential  equation  for  4  as  a  function 
of  p,  .  Therefore,  we  have  three  simultaneous  ordinary  differ¬ 
ential  equations  for 

A . (6.17) 

B. - (6.19) 

. . (6.27) 

and  one  algebraic  equation  for  A,,  namely  the  combination  of 
(6.20)  and  (6.27).  The  rest  of  the  coefficients  are  then 
determined  by 

(6.28)  At  «  3-3A .  -2 A,  -  £ 

(6.29)  Ai  *  f  -  a.  -  A,  -  A2 


(6.30) 

B, 

X  —  ^  Sm 

(6.31) 

Jx 

-  &  l  Be 

(6.32) 

s  *  J.6  Q 

(6.33) 

c, 

JJL  - 

^ fil  -  7.3718  * 

-  +  n  <,37i  -  0,*9SIX! 

fit 

4.  jo.  nt  _  10.06S1  +  wtS’X.' 
~PT~ 


(6.34) 


The  accuracy  of  these  polynomials  was  tested  for  the  Taylor 
case  where  the  exact  solution  is  known.  The  result  was  very 
encouraging;  for  Instance,  ft  was  given  accurately  by  the 
approximation  (6.27)  up  to  the  fourth  decimal  place.  This  does 
not  necessarily  mean  that  this  approximation  will  be  equally 
good  for  the  general  case,  but  it  shows  that  this  treatment 
gives  a  first  order  approximation  which  can  later  be  Improved 
b/  iteration. 
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6 . 5  Iteration  for  Improving  the  Approximate  Solutions .  The 
transformation  of  the  hydrodynamic  equations  (2.21)  to  (2.23) 
gives  the  solution  In  a  form  which  Is  relatively  Insensitive 
to  one  of  the  Independent  variables ,  namely^  .  This  suggests 
the  use  of  the  implicit  integrals  (4.8),  (4.11)  and  (4.15)  for 
iterations.  Once  good  approximations  Tor  p ,  J(  and  are  known, 
these  approximations  can  be  Improved  by  calculating  approximate 
>!($)  and  Hf)  for  (4.8)  and  £ ,  <f  ,  L  and  H  for  (4.10)  and  (4.11). 
Solving  these  for  f,  X  and  ^  would  yield  improved  approximations. 

J.  F.  Butler,  NOL,  has  derived  other  Integrals  which  can 
be  used  for  iterations  and  which  need  onlyjf  for  the  initial 
approximate  solution.  These  may  be  useful  because  the  average 

density  conditions  together  with  the  derivatives  x!  »  X,4 . 

lead  to  good  approximations  ofjf  without  much  numerical  effort. 

Solving  (2.22)  for  if  and  integrating  from  1  to  f  gives 

/ 


(6.35) 


The  reduced  pressure  can  be  obtained  from  (2.21) 


(6.36) 


-  6 


Numerical  integration  of  this  equation  yields  f($,  R>)  .  Using 
the  results  obtained  for  f  and  if  ,  (2.23)  can  be  used  to 
improve  jf  by  integrating  the  following  relation: 

(6.37)  #  -  ["■  * -A ¥}■ 


From  here  on  (4.5)  and  (6.36)  can  be  used  to  Improve  the 
functions  f  andjf  further. 
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CHAPTER  VII  THE  LOW  PRESSURE  RANGE 

When  />,  jpm  approaches  unity  it  becomes  sore  end  store 
difficult  to  evaluate  the  integral  #  because  J,  increases 
rapidly  with  decreasing  pt  andj*  .  For  values  of  close 

to  unity,  it  is  possible  to  derive  an  asysptotlc  relationship 
which  is  useful  in  solving  the  problem  in  the  range  of  low 
pressures.  The  first  step  in  doing  this  is  to  obtain  simple 
expressions  for  the  thermodynamic  data  of  water  which  permit  a 
representation  of  the  equation  for  <7,  in  dosed  form. 

7.1  Thermodynamic  Data  for  Water  at  Low  Pressures.  Por  low 
amplitudes,  the  Hugonlot  adiabatic  coincides  with  the  lsen- 
troplc;  therefore, 

n 

If  the  sound  velocity  changes  llnearlly  with  pressure, 

(7.2)  C  *  C0  (/  +  epj 

the  following  relations  can  be  derived: 


(7-3) 

?-/>• 

*  & 

/ 

/°° 

j  +  tp 

(7.4) 

»y 

=•  /  * 

tp 

(7.5) 

t 

=r  /  -h 

re  *  fa)  f> 

Strictly  speaking,  6  refers  to  the  change  of  the  sound  ve¬ 
locity  with  pressure  along  an  lsentroplc.  For  low  pressures, 
however.  It  is  permissible  to  use  values  of  €  which  are 
obtained  by  measuring  sound  velocities  as  functions  of  pressure 
at  constant  temperature,  namely  €  -  0.108.  This  holds  for 
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sea  water  as  well  as  fresh  water  at  all  temperatures  of 
interest.  For  low  pressures  (7*3)  to  (7*3)  yield  the  follow¬ 
ing  first-order  approximations: 


(7.6) 

“  A—/.  ' 

-  J  +OttX 

'**  p, 

p* 

(7.7) 

where 

/4^v- 

X  =r 

i 

*•  P» 

(p>  -Poj/p. 

«•  /  + 

^  (  for  fresh  water  at  20°  C 

<*.  -  -  2.376  f  f0i,  8ea  water  at  8°  C 


and  ql  ar  /  a  , 

Using  these  approximations ,  we  obtain  from  (3.9)  for  the 
dissipated  enthalpy  Increment 

(7.8)  ^  p.h  »-  />,*Z 

7.2  Behavior  of  the  Shockwave  Peak  Pressure.  As  the  sphere  of 
disturbance  Increases  in  size  the  shockwave  peak  pressure 
decreases  and  the  wave  finally  behaves  nearly  like  an  acoustic 
wave.  It  is  well  known  that  weak  shockwaves  never  attain  an 
exactly  acoustic  behavior  as  far  as  the  pressure  decay  with 
distance  is  concerned  i.e.,  the  pressure  does  not  exactly 
decrease  Inversely  proportional  to  distance.  For  instance, 
Kirkwood  and  Bethe  have  found  the  asymptotic  behavior  of  weak 
shockwaves  to  be  as  follows 

(7.9)  -A  «  $  a£j%jt 

ro  -  reference  length,  £  -  constant 


Ihla  gives 
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(7.10) 


*LL».r, 


~(>  *  til*. 


Our  expression  for  >  (3*7)  becomes  with  (7.6), 


(7.11) 


f>t  -~o 


where 


-J _ 

Z  *  /3  -  *  a  ~(  t  +  p) 

-  -(!  +  -t  ....  ) 

■  J  ,  ri  It *  jC~ 


comparison  of  these  two  equations  shows  that,  for  low  pressures, 
fi  must  slowly  (logarithmically)  approseh  unity.  This  means 
that  tends  to  zero,  when  pt  approaches  /o9. 

7.3  Expressions  for  9/  ,  x/  and  Y* '  at  Low  Pressures.  Using 
(7*6)  and  (7.7)#  we  obtain  the  following  limiting  equations  for 


H,  Jf,  end 

j  • 

Lmf  • 

^kVvw  .Hi. 

y,ltl 

(7.12) 

P*  po 

P.  -+P* 

~  -3 

'  -L-  9j_  -1 - 

'■  fc-  <*?„-&. 

/_=_OiJ. _ 

2  +f3  -Qr,xOt& 

■  Cl  - -3 JL  ...  ) 


+t»J 


(7.13) 


-&VW  L, 

P.  -/>.  ' 


■it/3 


H,  (t  -  2±x) 


Introduction  of  those  expressions  Into  the  equstlon  for  (p1 
yields 
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Ht. 


(7.14) 


s.  ,  , 

£  O  I  f*o 

'  (  (A  /J-"i 

<LX 


r-J 


-  X 


2-  L=A  -  -i-  2 —  /.  «*_  +  <*i  ) 
tt^X  J  +  /3  °j.  2*/i  1  *  TJ' 


Correspondingly,  we  obtain  for  %  '  : 


(7.15) 


J  jf,  !  l—£*  ( +  si.  j 

3,  -*/>.  '  l  Aa 


+  H,  ~hs_ 


°  z£o  ( 


’*  J 

P'M 


=  Oi-^xjj^~  ‘S/+2/  ~  ASlIA 

'  /  ->•/>.  J  2  +  ,i- a. 


—  ~^T7 

P.-/> + 


'=£ 
4 


^  A  IP,. 

2+/i 


__ 

P*~*  Pm 


£+■3 


+  ...  .  . 


for  the  derivative  of  the  reduced  density,  we  find 


(7 .16) 


A  -W>« 


*/ 


„  i  A~h  (y/ +  z  +  c,j\ 
4 -*■/>„  l  />•  '  J 


S  ^J-^VV 


7 A  The  Shockwave  Energy  and  its  Relation  to  the  Total  Energy. 
When  the  shockwave  has  propagated  to  large  distances,  the 
energy  distribution  within  the  sphere  of  disturbance  nay  be 
visualized  as  follows:  The  core  of  the  sphere  is  still  very 
hot*  Although  pressure  and  density  are  low,  the  region  near 
the  center  has  a  substantial  amount  of  internal  energy.  In 
the  case  of  underwater  explosions  this  energy  is  associated 
with  the  pulsations  of  the  bubble.  Near  the  surface  of  the 
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sphere,  l.e.,  near  the  shook  front,  one  finds  the  kinetic  and 
potential  energy  of  the  shockwave  aa  well  as  an  increment 
of  internal  energy  which  sterna  from  the  energy  dissipation 
at  the  shockfront.  Thus,  we  may  distinguish  three  energy  terms 
bubble  energy,  shockwave  energy  and  dissipated  energy.  In  the 
following,  we  deal  with  the  shockwave  energy. 

The  potential  energy  (reversible  Internal  energy)  of  the 
shockwave  is 


(7.17)  J/ 
where 

(7.18)  ^ 

For  low  pressures. 


K 

j C£if)  _  £/>„;  rLA*r 
c 

k 

-  fp  d  J  *  COHI  tj  •diP 


(7.19) 


If  we  Introduce  this  expression  into  the  Integral  (7.13)  severe 
terms  cancel  because  of  the  average  density  relation  (3.1), 

The  terms  left  are 

/ 


1/ 

A-V* 


+”•]£  T*t 

o 

/ 

K  P-fl,  f  n.  P-P, 

^  P*  J  fit 

O  , 


Afpf 


(7.20) 


flttf  «*  .4*1 
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where  (7.3)  has  been  uaed. 

Defining  a  reduced  shockwave  energy  9s  In  an  exactly  analo- 
goua  way  to  >)(  ,  (3.5)  and  ff  DUS  ,  (3.10),  namely: 

(7.21)  h,  -  —  ■  - - 


we  obtain 


(7.22) 


Pf  A 
*  -e/v 


r 

u  -  +  Jf  p's 


In  this  expression,  functions  for  »  X  and  V0  must  be  used 
which  refer  to  the  shockwave.  Suitable  functions  applicable 
to  the  low  pressure  range  are 


(7.23) 


X  "  f  „ 

f  -s 


This  yields 


**  *  X,' 
-  ?/ 


(7.24) 


Is  a  Pi 

/>•  -?• 


y  .*  ^  *9/ +3 


%'  +  X/  +  3 


or,  for  very  low  pressures,  neglecting  3  and  X,  ln  comparison 
with  2% 1  and  ^/; 


(7.25) 


*7  —  jL  d,  X  —  — 

A*  >V  <2  * 


It  Is  Interesting  to  note  that  this  result  can  be  also  obtained 
from  ^ 


(7.26) 


£pj  .  4n. rr3 J 


t-6o 


*  % 
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by  assuming  /ptejjfa  and  are  exponential  functions  , 
which  is  a  fairly  good  approximation  for  any  shockwave  when 
used  in  (7.26). 

Equation  (7 .25)  la  useful  only  if  an  inter-relationship 
between 


and 


/? « 

A*  * 


is  found.  For  this  purpose,  we  study  the  change  of  shockwave 
energy  with  the  distance  .  At  the  shock  front  mechanical 
energy  is  dissipated  due  to  the  irreversible  processes  associ¬ 
ated  with  the  discontinuous  rise  of  pressure,  velocity  and 
density,  as  discussed  in  Chapter  III. 

Hence 


(7.27) 


^cir,  sdt.r, 


47T  K  A 


Comparison  with  (3.11),  shows  immediately  that 

(7.28)  «  -  J  - — 


Combination  of  (7.26)  and  (3.7)  ylalda,  after  rearrangement, 

Zhfio  L ±_ 

*  ir=?o  ° *  ^ 

or  for  very  small  values  of  a 


(7.29) 


% 


(7.30) 


fiO 


Zhflt  2  +  /3 

p,x  /3».  -  ft 


Si.  x  LtA 
3  /3»."/3 
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0 


f 


Combination  with  (7.25)  gives  the  following  differential 
equation  for  : 

(7.3D  ~f  -T^to 

m  J  *Lt~\  X 

The  solution  of  (7*31)  is 


(7.32) 


foo 


[Si*  -  ^‘¥1  • 


is  the  integration  constant. 
Introduction  of  (7e32)  into  (7.25)  yields 


(7.33) 


/i  =  J  _  ^ 

/>-  -*l°o  '  J  tit 


with 

(7.34) 


A  rough  approximation  for  ft  la 


(7.35) 


fi  ~  -  JiL 

/>,  -»  />.  4t  ^ 


With  (7.11) >  we  obtain 

(7.36)  />,  —  />• 


-0  +  f 


/ 


This  expression  resembles  the  asymptotic  expression  of  Kirkwood 
and  Be the  (7.10).  For  the  actual  calculation,  approximation 
(7.35)  is  not  necessary.  It  is  almost  as  easy  to  use  (7.33). 
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7.5  The  Intermediate  Pressure  Range.  The  relationships  de¬ 
rived  In  Chapter  VII  are  suitable  for  large  values  of  ft,  ,  i.e., 
for  the  high  pressure  range.  Numerical  calculations  have  shown 
that  It  Is  difficult  to  extend  this  method  to  such  pressures , 
where  the  above  discussed  " low-pressure”  relations  become  valid. 
It  Is  possible  to  bridge  this  gap  by  means  of  the  "peak  ap¬ 
proximation" .  A  suitable  method  Is  described  In  reference  [h]. 
There,  two  ordinary  differential  equations  are  given  which 
relate  the  shockwave  peak  pressure,  time  factor  and  shape 
factor  with  distance.  The  time  factor  Is  given  by 


(7.37)  tm  ~ 


and  the  shape  factor  by 


(7.38) 


or, expressed  In  the  magnitudes  occurring  in  this  report: 

It- 


%  <b>  ( t  !&?/>,  *  Cr' 

-  &p(CK'*2lf  +ti)J 

In  these  equations  9*  and  c0  are  reference  values  which  cancel 
In  the  process  of  the  calculation.  Q0  may  be  chosen  as  an 
arbitrary  distance,  e.g.,  the  shock  radius  at  the  end  of  the 
high  pressure  range,  and  cQ  as  the  sound  velocity  of  the 


/  ,  / 


+ 


arfy'-r, 


/  2  1 
*t) 


(7-39) 


(A)  19* At 
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undisturbed  water.  Qio 7  and  Qsoy  are  derived  and  explained 
In  reference  [hj.  (7.37)and  (7.39)  as  well  as  the  peak 
pressure  calculated  a t  the  end  of  the  high  pressure  range, 
provide  the  Initial  conditions  for  the  peak-approximation 
method.  The  shape  factor  must  now  be  estimated  as  a  function 
of  distance.  This  is  not  difficult,  since  this  magnitude  Is 
close  to  unity  for  moderate  and  low  pressures.  A  function 
starting  at  the  calculated  values  in  the  high  pressure  region 
and  approaching  unity  Is  usually  sufficiently  accurate.  Inte¬ 
gration  of  the  differential  equations  for  peak  pressure  and 
time  factor  given  In  reference  [h]  makes  It  possible  to  obtain 
ft  for  such  low  values  of  pressure  that  (7*33)  becomes  valid. 
This  determines  the  Integration  constant  In  (7*34).  /3  can  be 

conveniently  expressed  In  the  terms  of  the  functions  ^  and 
which  occur  In  the  peak-approximation  theory: 


(7.40) 


/*h  A:-& 


l  8,  +  *£ 


This  provides  the  equations  for  the  calculation  of  the 
shocKwave  parameters,  beginning  at  the  extremely  high  pressures 
occurring  shortly  after  the  explosion  down  to  the  low  pressures 
of  a  nearly  acoustic  wave. 
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Figure  !• 


Determination  of  p  = 


c/Ai 


For  each  value  of  ft  several  values  of  @  are  assumed, 
sue*  as  0  *  i  1,  0,  -1,  -2,  etc.  For  each  of  these  values  of 
Q  n,  is  calculated  using  equation  1,3  and  In  is  plotted 
versus  In  .  The  correct  @  is  the  slope  of  the  solving 

curve  in  this  plot.  The  solving  curve  is  obtained  by  inter¬ 
polation  among  the  calculated  data.  Examples  of  this  data  for 
two  values  of  fj  are  shown  at  a  and  b  in  the  above  sketch. 

As  a  first  approximation,  a  curve  is  drawn  from  tie  known  initial 
value  of  at  ft  **  HfZ  through  the  regions  of  intersections  A 
and  B.  TV e  approximate  curve  is  then  adjusted  at  each  ft  , 
until  the  slope  of  the  curve  is  equal  to  the  interpolated  value 
oh  0  .  This  interpolation  is  readily  accomplished  with  the 
aid  of  a  crossplot  such  as  the  one  shown.  For  example,  the 
inclination  of  the  solving  curve  at  the  point  c  must  be  equal 
to  pt  read  from  the  P  scale  of  the  crossplot. 
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LIST  OP  SYMBOLS 

<X,  Qz  Abbreviations  defined  by  (7.6)  and  (7.7). 

A0)  A)  •  Coefficients  of  the  expansion  of  V' 

B0f  8/  .Coefficients  of  the  expansion  of  X 
C  Sound  velocity 

c0  Sound  velocity  of  the  undisturbed  medium 

C0j  ....  Coefficients  of  the  expansion  of 

t  2.71828 _ 

E  Internal  energy 

^  Abbreviation  defined  by  (4.8) 

3  Abbreviation  defined  by  (4.13a) 

6,  Decay  factor  of  shock  front  density  (2.20) 

h  Dissipated  enthalpy  Increment 

AH  Enthalpy  increment 

H  Abbreviation  defined  by  (4.9) 

//,  Decay  factor  of  shockwave  peak  pressure  (2.18) 

J  Reduced  internal  energy  (2.9' 

Abbreviations  defined  by  (4.13) 

^  Abbreviation  defined  by  (4.9) 

L,  Decay  factor  of  partial  velocity  behind  shock 
front  (2.19) 

m  m  ) 

Abbreviations  defined  by  (5*9) 

Excess  pressure 
/p,  Shockwave  peak  pressure 

Po  Absolute  pressure  in  the  undisturbed  medium 
<J  Abbreviation  defined  by  (4.8) 

Q  Energy  yield 

Dissipated  energy 
Qg  Bubble  energy 
Shockwave  energy 
r  Radius 
5  Entropy 
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t  Time 

7  Temperature 

T*  Lowest  temperature  of  the  system 
ac  Particle  velocity 

A*  Particle  velocity  right  behind  shock  front 
U  Propagation  velocity  of  shock  front 
v  Specific  volume  •  '/A 

X  m  (ft -/>*)/ ft* 

j 

J  Adiabatic  exponent  as  defined  by  (2.4) 

True  adiabatic  exponent;  see  (2.5) 

^  Adiabatic  exponent  right  behind  the  shock  front 

Adiabatic  exponent  at  the  center  of  the 
v  disturbance 

Ranklne-Hugonlot  exponent 

4,  ii  Roots  of  denominator  in  (4.11) 

5  Coefficient  of  change  of  sound  velocity  with 
pressure 

Reduced  energy 

?/  Total  reduced  energy  of  the  sphere  of  disturbance 
^4,  Reduced  dissipated  energy 
Reduced  shockwave  energy 

A  -  ~  AVA 

?  Reduced  radius 
7T  3.14159.... 
p  Density 

f>,  Density  right  behind  shock  front 

/?,r  Density  rlxht  behind  shock  front  at  the  initial 
r  condition  (Taylor  Solution) 

Po  Density  of  the  undisturbed  medium 

^  Reduced  particle  velocity 

X  Reduced  density 
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Reduced  pressure 

Reduced  pressure  at  the  center  of  the 
disturbance 
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